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Abstract. Nonanalyticities of thermodynamic functions are studied by adopting an approach based 
on stationary points of the potential energy. For finite systems, each stationary point is found to cause 
a nonanalyticity in the microcanonical entropy, and the functional form of this nonanalytic term is 
derived explicitly. With increasing system size, the order of the nonanalytic term grows, leading to 
an increasing differentiability of the entropy. It is found that only "asymptotically flat" stationary 
points may cause a nonanalyticity that survives in the thermodynamic limit, and this property is 
used to derive an analytic criterion establishing the existence or absence of phase transitions. We 
sketch how this result can be employed to analytically compute transition energies of classical spin 
models. 
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A phase transition is an abrupt change of the macroscopic properties of a many- 
particle system under variation of a control parameter. An approach commonly used 
for the theoretical description of phase transitions is the investigation of the analyticity 
properties of thermodynamic functions like the canonical free energy of enthalpy. It is 
long known that nonanalytic behavior in a canonical or grandcanonical thermodynamic 
function can occur only in the thermodynamic limit in which the number of degrees of 
freedom N of the system goes to infinity [1]. Recently, however, it was observed that the 
microcanonical entropy, or Boltzmann entropy, sn of a finite system is not necessarily 
real-analytic, i. e. not necessarily infinitely many times differentiable. More specifically, 
it was observed in [2] that stationary points of the potential energy function V^(q) of a 
classical many-body system with continuous variables q = (qi , . . . , q^) give rise to these 
nonanalyticities. Generically, with increasing N, the nonanalyticities appear in higher 
and higher derivatives of the microcanonical entropy spj. Despite this "smoothening" of 
the entropy when approaching the thermodynamic limit, it was shown recently that a 
finite-system nonanalyticity of sn can survive the thermodynamic limit if the Hessian 
determinant of Vn, evaluated along a suitable sequence of stationary points for different 
system sizes N, goes to zero in a suitable way. This observation was used in [3, 4] to 
derive an analytic criterion, local in microscopic configuration space, on the basis of 
which the existence of phase transitions can be analyzed and, in some cases, an exact 
analytical expression for the phase transition energy can be derived. 

In this article, the functional form of the nonanalyticities of the finite-system entropy 
sn is reviewed. Based on this result, we then sketch the criterion which relates the 
occurrence of a phase transition to the vanishing of the Hessian determinant evaluated 
along a sequence of stationary points. Finally, a simple strategy is discussed which 




FIGURE 1. As a simple example, consider a two-dimensional configurational space R 2 on which a 
potential energy function ^(^1,^2) = <7i +I2 ^ s defined (left). The corresponding density of states as 
defined in (2) has a discontinuity at v = vq (right). 



permits to construct sequences of stationary points for systems of arbitrary size N. 



NONANALYTICITIES OF THE FINITE-SYSTEM ENTROPY 

Canonical and grandcanonical ensembles are the ones most frequently used in statisti- 
cal mechanics applications. The corresponding thermodynamic potentials, i.e. the free 
energy and the grandcanonical potential, are known to be analytic functions for all fi- 
nite systems sizes, and many physicists have become so used to this property that they 
expected all finite-system thermodynamic functions to be smooth. It is, however, fairly 
easy to construct counterexamples to this false expectation. The thermodynamic func- 
tion we consider here is the configurational microcanonical entropy 1 

s N (v) = lnQ N (v)/N, (1) 

where 

a N (v) = J dq8(V N (q)-Nv) (2) 

is the configurational density of states. The integration in (2) is over configuration space, 
and v denotes the potential energy per degree of freedom. For a potential energy function 
V2(qi,q2) = <?i +<?2' tne configurational density of states ^(v) is easily shown to have 
a discontinuity at v = vo (see figure 1). For a slightly less trivial example of a potential 
with a proper saddle point, see figure 1 of [6]. 

In both examples, the nonanalyticities of Q.n(v) occur precisely at the values of the 
potential energy per degree of freedom v s = Vjv(q s )/N which correspond to stationary 
points of the potential, i.e. points q s where dV/v(g s ) = 0. This observation remains valid 
in higher dimensional configuration spaces, and independent of whether the stationary 
point is a minimum, a maximum, or a saddle point. In the following, in order to obtain 
a general result characterizing the nonanalytic behavior of Q.^ induced by a stationary 
point, we require the potential Vn to be a Morse function, i.e. to have a non-vanishing 
determinant of the Hessian at all stationary points of Vn- One may argue that this 



Nonanalyticities of the "usual" microcanonical entropy are related to nonanalyticities of its configura- 
tional counterpart in a rather straightforward way. See [5] for details. 



is an insignificant restriction, since Morse functions form an open dense subset of the 
space of smooth functions [7] and are therefore generic. This means that, if the potential 
Vn we are interested in is not a Morse function, we can transform it into one by adding an 
arbitrarily small perturbation. An important consequence of the Morse property is that 
all stationary points of such a function are isolated which allows us to study the effect of 
a single stationary point on the configurational density of states (2). Such an asymptotic 
analysis has been reported in [8, 4]: 

Theorem 1 Let Vm '■ G — > R be a Morse function with a single stationary point q s of 
index k in an open region G C H N . Without loss of generality, we assume V^(q s ) = 0. 
Then there exists a polynomial P of degree less than N/2 such that at v = the 
configurational density of states (2) can be written in the form 

a N (v) = P(v) + - hN ' k(v) + o{v N / 2 - £ ) (3) 

for any £ > 0. Here is the Heaviside step function, o denotes Landau's little-o symbol 
for asymptotic negligibility, and 

{(_l)V2 v (A'-2)/2 ^ fork even, 

(_l}(k+\)/2 v {N-2)/2 % -\ Xn \ v \ for N even, k odd, (4) 
( _ 1 ) {N-k)/2 ( _ v ) (JV-2)/2@( _ v ) for ^ k oM 

For a proof of this result see [4]. In short, we see from theorem 1 that, at a nonanalyticity 
of £2;v(v) induced by the presence of a stationary point of Vm, the configurational 
density of states is [(N — 3) /2J -times differentiable at the corresponding value v s of 
the potential energy. Hence, when increasing the number /V of degrees of freedom, the 
non-differentiability occurs in higher and higher derivatives of Q#(v) [or Syv(v)], and 
one might naively expect the nonanalyticity to disappear in the thermodynamic limit. 



FLAT STATIONARY POINTS AND PHASE TRANSITIONS 

The result of the previous section, and in particular the unbounded growth of the differ- 
entiability with increasing N, does not seem to suggest any connection between station- 
ary points of the potential energy function and phase transitions in the thermodynamic 
limit N — > oo. There is, however, a result by Franzosi and Pettini which indeed estab- 
lishes such a relation. The following sloppy reformulation of the result will be sufficient 
for our purposes. 

Sloppy theorem 1 Let be the potential of a system with N degrees of freedom 
and short-range interactions. If some interval [a,b] of potential energies per degree 
of freedom remains, for any large enough N, free of stationary values of Vn, then the 
configurational entropy s(v) = limA?^ 00 5^(v) does not show a phase transition in this 
interval. 

Note that a precise formulation of this theorem requires further technical conditions on 
the potential Vn (see [9, 10] for details). 



Indeed, the finite-system result of theorem 1 can be helpful towards an understanding 
of how nonanalyticities of the entropy s^ can give rise to a nonanalyticity in the thermo- 
dynamic limit: From equation (3), we observe that the nonanalytic term h N ^ comes with 
a prefactor of 1 / y/\det [J^ v (q s )]\. Although, with increasing N, the differentiability of 
hffk grows unboundedly, this "smoothing" may be counterbalanced by a vanishing (with 
increasing N) Hessian determinant. This intuition has been made rigorous in [3, 4], but 
here we will give only a sloppy reformulation capturing the essence of this result. 

Sloppy theorem 2 The nonanalytic contributions of the stationary points of V/v to the 
configurational entropy cannot induce a phase transition at a potential energy per 
particle v t if in a neighborhood ofv t , 

1. the number of critical points is bounded by exp(CAf) with some C > 0, and 

2. the stationary points do not become "asymptotically flat" in the thermodynamic 

limit, i.e. lim_/v^ 00 |det e ^(^ s ) \ l ^ N is bounded away from zero for any sequence of 
stationary points q s lying in the vicinity ofv t . 

For a precise formulation and a proof of this result see [4]. In short, this result classifies 
a subset of all stationary points of V/v as harmless as what regards phase transitions 
and leaves only the asymptotically flat ones as candidates capable of causing a phase 
transition. 



SPECIAL SEQUENCES OF STATIONARY POINTS 

Importantly for the application of sloppy theorem 2, knowledge of a suitably chosen 
subset of the stationary points of V/v may be sufficient: If one manages to find some 
sequence of stationary points such that, along this sequence, 

lim\detJl? v (q s )\ l/N = 0, (5) 

the corresponding limiting value v t = lim/v^co V/v (q s ) /N is a good candidate for the exact 
value of the phase transition potential energy. This idea was first employed by Nardini 
and Casetti in [1 1], where suitably constructed sequences of stationary points were used 
to single out the phase transition of a model of gravitating masses and determine its 
critical energy. 

To illustrate how special sequences of stationary points can be constructed, we con- 
sider a one-dimensional XY model with periodic boundary conditions, characterized by 
the potential energy function 

v»(,)^T '-^- w) 

i=\ j=i J 

where qi e \—K,%) are angular variables, and a is some nonnegative exponent. For 
a G [1,2], this model is known to show a phase transition from a ferromagnetically 
ordered to a paramagnetic phase, but no exact thermodynamic solution is known. 




FIGURE 2. Sketch of stationary points of Vn for N = 8, where is the angle between the arrow and 
the dashed axis. Left: Stationary points where all qi e {0,7c}. Right: All differences qk — qk-\ between 
neighboring angles are equal, with differences chosen such that qo = q^, in compliance with the periodic 
boundary conditions. 



Stationary points of the potential energy (6) have to satisfy the set of equations 

Q = dV N (q) = (yv ^, )/2 sulfas - q k+j ) + sin(q k - q k _j) 
dqk jTi J a 

for k = 1, ... ,N. To get rid of the trivial global rotational invariance of (6), we fix = 
and eliminate the equation with k = N in (7). The thermodynamics of this reduced model 
is identical to that of the full one, as the contribution of one degree of freedom to the 
partition function is negligible in the thermodynamic limit. 

There are two particularly simple classes of solutions of (7), similar in spirit to 
those constructed in [11] for a one-dimensional model of gravitating masses: First, any 
combination of qi e {0, n} for i = 1, . . . ,N — 1 will make the sine functions in (7) vanish. 

A second class of solutions is given by q$ = 2nmn/N for m,n e {l,...,Af}. These 
solutions have equal angles between neighboring spins. As a result, sin(<7^ — qk+j) = 
sm(q k _j : — q k ), and therefore each of the summands in (7) vanishes separately. Both 
classes of solutions are sketched in figure 2. To employ these classes of stationary points 
along the lines of sloppy theorem 2, one needs to evaluate the Hessian determinant of 
(6) at the stationary points. This is work in progress and will be reported elsewhere. 
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